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Abstract. We present a detailed study of Higgs boson production in association with a single top quark 
at the LHC, at next-to-leading order accuracy in QCD. We consider total and differential cross sections, 
at the parton level as well as by matching short distance events to parton showers, for both t-channel and 
s-channel production. We provide predictions relevant for the LHC at 13 TeV together with a thorough 
evaluation of the residual uncertainties coming from scale variation, parton distributions, strong coupling 
constant and heavy quark masses. In addition, for t-channel production, we compare results as obtained 
in the 4-flavour and 5-flavour schemes, pinning down the most relevant differences between them. Finally, 
we study the sensitivity to a non-standard-model relative phase between the Higgs couplings to the top 
quark and to the weak bosons. 


1 Introduction 

The first Run of the LHC has already collected compelling 
evidence that the scalar particle observed at 125 GeV is the 
one predicted by the Brout-Englert-Higgs symmetry breaking 
mechanism of SU{2)l x U{1)y as implemented in the 
Standard Model (SM) [^. In such minimal case, the strengths 
of the Higgs boson couplings to the elementary particles, in¬ 
cluding the Higgs boson itself, are uniquely determined by their 
masses. While somewhat limited and subject to additional ad 
hoc assumptions, the first measurements of the Higgs couplings 
to fermions and vector bosons agree well with the SM predic¬ 
tions [^[^. 

Such general agreement with the SM expectations and the 
absence of any evidence (from the LHC itself) of the existence 
of new states at the TeV scale, motivate a thorough study of 
the Higgs boson interactions at the Run H. In addition to the 
coupling strength determinations conducted so far, the Lorentz 
structure of the vertices as well as the possible existence of a 
relative phase among the couplings need to be fully assessed. 
In order to gather the necessary information, the widest possi¬ 
ble campaign of measurements has to be undertaken, including 
different production and decay modes of the Higgs boson. In 
addition, given the limited discriminating power of single chan¬ 
nels, a global combination of the relevant measurements will be 
necessary. To achieve this goal at the LHC one needs to adopt a 
complete and consistent theoretical framework, able to encom¬ 
pass interactions that go beyond the SM (and possibly to or¬ 
ganise them in terms of an ordering principle), and that allows 
the systematic inclusion of higher-order corrections, both QCD 
and electroweak (EW). This latter point is a conditio-sine-qua- 
non at the LHC, in order to control total rates and differential 
distributions and to estimate the residual uncertainties. Such a 


theoretical framework exists and amounts to “simply” extend 
the dimension-4 SM Lagrangian to all operators of higher di¬ 
mensions (up to dimension -6 in this first instance) consistent 
with the unbroken SM symmetries SU{3)c x SU{2)l x U{1)y', 
i.e. to consider the SM as an effective field theory valid up to 
a scale A [^[^. 

This work fits in the above general strategy and focuses 
on Higgs production in association with a single top quark. 
As in single top production, at the leading order (LO) in QCD 
one can organise the production mechanisms into three groups, 
based on the virtuality of the W boson: t-channel production 
(fig- 0 - s-channel production (fig. 0 , and associated produc¬ 
tion with an on-shell W boson. While characterised by a rather 
small cross section with respect to the main single Higgs pro¬ 
duction channels (gluon-gluon fusion, vector boson fusion and 
associated production, and tiH), Higgs and single-top associ¬ 
ated production features unique aspects that make this pro¬ 
cess particularly interesting for Higgs characterisation 00. 
Notably, it is among the very few processes relevant for LHC 
phenomenology (together with 77 —> 77 and gg —>■ ZH) to be 
sensitive to the relative size and phase of the coupling of the 
Higgs boson to the top quark and to the weak bosons. For t- 
channel and W-boson associated production, diagrams where 
the Higgs couples to the top quark interfere destructively with 
those where the Higgs couples to the W boson (due to the 
unitarity of the weak interactions in the SM), making cross 
sections and distributions extremely sensitive to departures of 
the Higgs couplings from the SM predictions |10| . 

The aim of the hrst part of this work is to provide accurate 
SM predictions including QCD corrections at next-to-leading 
order (NLO) for t- and s-channel Higgs production in associ¬ 
ation with a single top quark, as well as reliable estimates for 
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Fig. 2. LO Feynman diagrams for s-channel tH production. 

the residual uncertainties in rates and distributions. Particular 
attention is devoted to the uncertainty related to the different 
flavour schemes that can be adopted to compute the dominant 
t-channel production mode. The corresponding SM predictions 
are the necessary theoretical input to possibly assess the exis¬ 
tence of deviations due to new physics (be it resonant or not); 
to this aim, the study of the uncertainties in total rates as well 
as in differential distributions becomes of primary importance. 

We then consider how accurately and precisely the effects 
of the (only) dimension-6 operator that modifies the value and 
the phase of the top quark Yukawa coupling can be predicted, 
again at the total as well as at the differential level. This in¬ 
formation is useful to assess the reach of the LHC to constrain 
the relevance of this dimension-6 operator {i.e. to bound the 
complex coefficient in front) and, if deviations from the SM are 
detected, to quantify them. 

The paper is organised as follows. In sect, [^we introduce 
the main features of the Higgs and top quark associated pro¬ 
duction. In sect. 1^ we focus on the t-channel production mode, 
with a special attention to the issues connected to the 4-flavour 
(4F) and 5-flavour (5F) schemes. We describe the settings of 
the calculation, present results in the SM for total rates up to 
NLO in QCD and their uncertainties, and finally show relevant 
differential distributions at NLO matched to a parton shower. 
In sect. we shortly consider the s-channel production mecha¬ 
nism, which has a much smaller impact on Higgs phenomenol¬ 
ogy in the SM. We evaluate the total cross sections and its 
uncertainties in the SM and show some representative distri¬ 
butions in comparison with the corresponding t-channel ones. 
In sect. we study the impact of an anomalous, CP-violating 
top quark Yukawa interaction on t-channel production, both 
at the total and differential cross section level. We summarise 
our findings in sect. 


2 Main features 


organise the various production mechanisms into three groups, 
based on the virtuality of the W boson: f-channel production 
features a space-like W, s-channel production a time-like W, 
and VF-associated production an on-shell W boson. One has 
to bear in mind that while this classification is certainly use¬ 
ful, it is not physical, being an approximation that holds only 
at LO and in the 5-flavour scheme. At higher orders in QCD, 
or using a different flavour scheme to define the processes, the 
separation becomes increasingly fuzzy, as it will be clarified at 
the end of this section. 

As in single top production in the SM, tH production is 
always mediated by a tWh vertex and therefore it entails the 
presence of a bottom quark either in the initial (t-channel and 
W-associated) or in the final state (s-channel). In the case of 
initial-state bottom quarks, two different approaches, the so- 
called 4F and 5F schemes, can be followed to perform pertur¬ 
bative calculations. 

In the 4F scheme one assumes that the typical scale of 
the hard process Q is not significantly higher than bottom 
quark mass, which in turn is considerably heavier than Aqcd, 
Q ^ ^ Aqcd- Technically, one constructs an effective the¬ 

ory of QCD with only four light flavours, where heavier quarks 
(bottom and top), being massive, do not contribute to the 
initial-state proton wave-function (in terms of parton distribu¬ 
tion functions (PDFs)), nor to the running of the strong cou¬ 
pling, and they appear only as final-state particles. In so doing, 
mass effects in the kinematics of heavy-quark production are 
correctly taken into account already at the lowest order in per¬ 
turbation theory. In addition, the matching to parton-shower 
programs is straightforward, the heavy-quark mass acting as an 
infrared cutoff for inclusive observables. However, limitations 
might arise when Q 3> mb and one probes kinematic configu¬ 
rations which are dominated by almost collinear g ^ bb split¬ 
tings: in this case the accuracy of predictions can be spoiled by 
large logarithms log{Q^/mf) appearing at all orders in pertur¬ 
bative QCD. Were this the case, such large logarithms would 
harm the behaviour of a fixed order expansion in as- 

This issue can be addressed in the 5F scheme (and improve¬ 
ments thereof), whose aim is to reorganise the perturbative ex¬ 
pansion by resumming such logarithms via the DGLAP equa¬ 
tions. One starts by assuming Q ';$> mb and defines a scheme 
where power corrections of order mbjQ appear at higher orders 
in the as expansion. In practice, one sets the bottom mass to 
zero and includes bottom quarks in the initial state as proton 
constituentsj^In so doing, towers of logarithms associated with 



In this section we introduce the main features of Higgs produc¬ 
tion in association with a single top quark. As already men¬ 
tioned in the introduction, at LO in QCD one can effectively 


^ The bottom mass can be reinstated explicitly at higher- 
orders by systematically including it in diagrams that do not 
feature bottom quarks in the initial state, the so-called S- 
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the initial-state g ^ bb splitting are resummed to all orders in 
perturbation theory by evolving the perturbative bottom quark 
PDF via the DGLAP equations. 

Computations in the 5F scheme are typically much simpler 
than the corresponding 4F ones, because of the lesser final-state 
multiplicity and the simpler phase space. This is for example 
the reason why single-top production is known at NNLO in the 
5F scheme while only at NLO in the 4F [^. For a system¬ 
atic investigation of the sources of differences between the 4F 
and 5F schemes in single fe-quark and double fe-quark induced 
processes we refer the reader to [14| and [15| , respectively. In 
short, the 4F and 5F schemes differ in what kind of terms 
are pushed into the missing higher-order corrections. There¬ 
fore, as the accuracy of the predictions for a given observable 
increases, milder differences should be expected between the 
schemes. This provides a strong motivation to go at least to 
NLO accuracy in the computation of the t-channel cross sec¬ 
tion, in order to reduce the flavour-scheme dependence of the 
predictions and thus the overall theoretical uncertainty. The 
final accuracy, however, will depend on the specific observable 
considered, whose perturbative accuracy can be different in the 
two schemes. 

In the case of (Higgs and) single top production at hadron 
colliders, the 5F scheme has also the operational advantage 
that allows an easy separation of the various production mecha¬ 
nisms into the three groups mentioned above. In the 5F scheme 
the t-channel, s-channel and IT-associated production are in¬ 
dependent up to NLO and start to interfere only at NNLO, 
and the kF-associated production interferes with tiH starting 
from NLO. In the 4F, on the other hand, the t-channel at 
NLO can interfere with the s-channel (at NNLO) and with 
kF-associated production (if the kF decays hadronically), and 
the kF-associated production also interferes with tiH already 
at the tree level. While the former interferences are very small 
and can be safely neglected if the aim is to evaluate the dom¬ 
inant t-channel cross section, the interference of kF-associated 
production with tiH turns out instead to be quite large. The 
on-shell kF associated production therefore needs a dedicated 
study that we defer to a separate work. 


3 t-channel production 

In this section we present the SM predictions for f-channel 
Higgs plus single top production at the LHC (see fig. [^, at 
NLO accuracy in QCD. We first describe the technical setup 
we have used for NLO simulations, the input parameters as 
well as the various sources of theoretical uncertainties. We 
then show results for the inclusive cross section at the LHC 
with y/s = 13 TeV, discussing how to combine the theoretical 
uncertainties, and finally present NLO distributions matched 
to parton shower. 


3.1 NLO simulations, parameters and uncertainties 

In this work, we employ the MadGraph5_aMC@NLO frame¬ 
work [^, which allows to compute both inclusive cross sec¬ 
tions and differential distributions matched to parton-shower 

ACOT scheme . In this work we adopt a “pure” 5F scheme 
where mb — 0 throughout. 


programs, up to NLO accuracy in QCD, in a fully automatic 
way 17 -20 once the relevant Feynman rules and UV/i ?2 coun¬ 


terterms for a given theory are provided in the form of a UFO 
model 21 ■ 23 . While these extra Feynman rules are available 
in MadGraph5_aMC@NLO by default for the SM, non-SM 
interactions that will be considered later in sect. [5] are encoded 
in the HC_NLO_XO model |24f 26 , publicly available online in 
the FeynRules repository [27] . 

In MadGraph5_aMC@NLO the code and events for t- 
channel tH production at hadron colliders in the 4F scheme can 
be automatically generated by issuing the following commands: 


(> import model loop_sm) 

> generate p p > h t b~ j $$ w+ w- [QCD] 

> add process p p > h t~ b j $$ w+ w- [QCD] 

> output 

> launch 


while the corresponding commands in the 5F scheme are: 

> import model loop_sm-no_b_mass 

> define p = p b b~ 

> define j = p 

> generate p p > h t j $$ w+ w- [QCD] 

> add process p p > h t~ j $$ w+ w- [QCD] 

> output 

> launch 


Note that the $$ w+ w- syntax removes s-channel tH diagrams 
as well as real-correction diagrams where an on-shell IF decays 
to two light quarks, which belong to kF-associated produc¬ 
tion. The top quark decays are subsequently performed start¬ 
ing from the event file (in the Les Houches format 28 ) by 
MadSpin [29] , following a procedure 30 that keeps spin cor¬ 
relations. 

In the numerical calculation, the mass of the Higgs boson is 
set to mrr = 125.0 GeV, while the mass of the top quark is set 
to mt = 173.3 GeV. We renormalise the top quark Yukawa cou¬ 
pling on-shell, setting it to ytl\/2 = mt/v, where v ~ 246 GeV 
is the EW vacuum expectation value. 

The on-shell mass of the bottom quark is set to 


mt, = 4.75 ±0.25 GeV, 


( 1 ) 


where we take the uncertainty to be of 0{Aqci>), accordingly 
to the prescription in ref. 31 . On the other hand, we set the 


bottom quark Yukawa coupling to zero, because effects related 
to the Hbb interactions are negligible for this process. We re¬ 
mind that in the 4F scheme the value of mb enters the hard- 
scattering matrix element and the final-state phase space, while 
in the 5F scheme it affects only the parton luminosity. 

PDFs are evaluated by using three global fits: NNPDF2.3 132] 


MSTW2008 33 and CTIO 34 , through the LHAPDF inter¬ 


face [35]. PDF uncertainties are computed for each PDF set. 


following the recipes summarised in 36 . A comparison among 


these three global fits allows to estimate the PDF systematic 
uncertainties related to the technical details of the fitting pro¬ 
cedure employed by each group. We note that the above three 
PDF collaborations provide NLO PDF sets both in the 4F 
and 5F schemes, while only MSTW gives LO PDFs in both 
the schemes. 

The reference value for the strong coupling constant we 
employ here is 


= 0.1190 ±0.0012, 


( 2 ) 
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Fig. 3. Scale dependence of the total cross sections for the pp —> tHq + tHq production at the 13-TeV LHC, where the 4F 
(blue) and 5F (red) schemes are compared. LO (dashed) and NLO (solid) predictions with MSTW2008 LO/NLO PDFs are 
presented for pjj = pp = p, with a static (left figure) and a dynamic (right figure) scale choice. Two off-diagonal profiles of the 
scale dependence at NLO are also shown, for (pr? = V^p, pp = pjy/^') and for (prr = p/%/2, pp = %/2p). The black arrows 
visualise the envelope of the combined scale and flavour-scheme uncertainty defined in eq. 


where the uncertainty is taken accordingly to the PDF4LHC 
recommendation [36[|37| , and the central value is chosen such 
that our 68% confidence interval encompasses the current PDG 
world average and the best as{mz) estimates obtained by 
each of the three PDF global fits |39}|4l| . We remark that the 
value in eq. ([^ is consistent with the 5F description. Since the 
difference between 4F and 5F in the running is limited to 
scales above mt, eq. ([^ can be translated into the following 
condition on as(mi,) (running as at 2-loop accuracy) 

Qr^°^(m6) = 0.2189 ± 0.0042 , (3) 


which is now flavour-scheme independent. 

CTIO does not provide PDF sets to compute mt uncertain¬ 
ties in the 5F scheme and PDF uncertainties in the 4F scheme; 
both CTIO and MSTW2008 do not provide 4F PDF sets with 
different as(mz) values. Thus, it is possible to address all the 
various sources of uncertainty in both schemes only when using 
NNPDF2.3 parton distributions, while MSTW2008 and CTIO 
uncertainty bands can be sometimes underestimated (though 
just slightly, as we will see later in sect. |3.2| ). 

For matching short-distance events to parton shower we use 
the MC@NLO method with Pythia8 [^, while HER- 
WIC6 [43| has been used for a few comparisons. We recall that 
matching to Pythia6 44 (virtuality-ordered, or pr-ordered 
for processes with no final-state radiation) and HERWICH—h 
1451 are also available inside MadCraph5_aMC@NLO. Jets 


are reconstructed by means of the anti-fcr algorithm 46 


plemented in FastJet [47| , with distance parameter R = 0.4, 
and required to have 


jet satisfies 

Prijb) > 30 CeV, \v{jb)\ < 2.5. (5) 

We assume 100% b-tagging efficiency in this work. 


3.2 Total rates 

In this section we present the total cross section for t-channel 
production of a Higgs boson together with a single top quark 
(or antiquark), at NLO in QCD. The main sources of theoret¬ 
ical uncertainty that we address here are: 

— renormalisation and factorisation scale dependence, 

— 4F and 5F scheme dependence, 

— PDF uncertainty, 

— asijnz) uncertainty, 

— rtib uncertainty. 

At the end of this section we will also briefly comment on the 
impact of the bottom quark Yukawa coupling and of the de¬ 
pendence of the results on the Higgs and the top quark masses. 

We start by showing in fig. the renormalisation and fac¬ 
torisation scale dependence of the LO and NLO total cross 
sections, both in the 4F and 5F schemes. We compute cross 
sections with two different scale choices, and vary pp = pp = p 
around a central scale po which is chosen as 

pI = {tuh + ( 6 ) 


PtU) > 30 GeV , |p(j)| < 4.5 . 


(4) 


for the static scale choice (left figure), and 


A jet is identified as b-jet if a b-hadron (or b-quark for fixed- Mo — Hp/G — E mT{i)/6 (7) 

order calculations) is found among its constituents, and if the i^H,t,b 
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scheme 

ULO [fb] UNLO [fb] K 

4F (ho) 

5F (m?) 

63.46(8) 69.43(7) 1-09 

60.66(6) 73.45(8) 1-21 

4F ini) 
5F ini) 

64.31(8) 71.29(10) 1.11 

58.83(5) 71.54(7) 1.22 


Table 1. LO and NLO cross sections and corresponding K 
factors for t-channel tH production at the 13-TeV LHC in the 
4F and 5F schemes. MSTW2008 PDFs have been used. The 
integration error in the last digit (s) and the scale dependence 
by a factor 2 around the static and dynamic scale choices in 
eqs. § and Q are also reported. 


for the event-by-event dynamic choice (right figure), where 
rriT = \Jm? -f is the transverse mass of a particle. 

We find a pattern similar to the case of the single top pro¬ 
duction (see fig. 3 in [^). At LO the scale dependence in the 
4F scheme is stronger than in the 5F, simply because the 4F 
calculation starts already at order as. As expected, predictions 
at NLO are much more stable under the scale variation than 
at LO. We find that the 4F and 5F predictions are in better 
agreement if /r is chosen to be roughly a factor 4 ( 6 ) smaller 
than the typical hard scale of the process thh + rtit (Ht) for 
the static (dynamic) scale choice. This is a known and gen¬ 
eral feature of 6 -initiated processes at hadron colliders . At 
such reduced scales the 4F and 5F predictions are typically in 
good agreement, and this is indeed what we observe taking the 
reference scale choice po as in eqs. § and Q. Table shows 
the corresponding values of the LO and NLO cross sections 
in fig. i where the uncertainty from missing higher orders is 
estimated varying the scale /r by a factor 2 around fj,o- 

In fig. 11 we also plot two off-diagonal (p_R 7 ^ fip) slices of 
the NLO cross section surface in the plane {fiu, fip), shifted 
by a factor in the direction orthogonal to the diagonal. 
The effects of off-diagonal scale choices are more pronounced 
in the 4F scheme than in the 5F, even though in general they 
are quite modest, except at very low scales, i.e. comparable to 
mb- We conclude that, for our choice of fio, the diagonal hr = 
Hf IS sufficiently representative of the scale dependence of the 
total cross section, when the scale is varied by the usual factor 
two. We also observe that the scale value which minimises the 
flavour-scheme dependence is rather stable under shifts away 
from the diagonal. 

We note that the scale dependence pattern is strongly cor¬ 
related to the flavour scheme employed. Therefore, after we es¬ 
timate the scale dependence of both 4F and 5F results (varying 
the scale hf = fJ^R = /r by a factor 2 around fio), we define a 
combined scale and flavour-scheme uncertainty band by tak¬ 
ing the envelope of the extremal points (shown by the black 
arrows in fig. , and the best prediction for the cross section 
as the central point of this envelope. The total cross section at 
NLO and its combined scale plus flavour-scheme uncertainty 
are defined by 


o-NLO = (u-’''-I-n )/2, 5 ^+fs = (cr’^ - cr )/2 , (8) 


where 

(T+= max {crN'Lo(M) , o-NLo(p)} . (9) 

Me[/io/2,2Mo] 

a~ = min > crNLo(M)} ■ ( 10 ) 

Mebo/2,2Mo] 

Now we turn to the PDF, as{mz) and mb uncertainties. 
In principle these three uncertainties can be correlated. How¬ 
ever, the correlations are very small and can be often neglected 
in combinations. For example, using NNPDF, we have explic¬ 
itly checked that the combined PDF-fOs uncertainty computed 
with full correlations differs from the one without correlations 
by 0.1% at most. In the 4F scheme mb is independent of PDF 
and Os, while we confirmed that the uncertainty correlation 
between PDF and mb in the 5F scheme is well below the per¬ 
cent level. Moreover, the correlation between a^ and mb is tiny 
and can be neglected . We note that neglecting correlations 
allows us to compare PDF uncertainty bands at a common Oa 
value, once central predictions (computed with this common 
tts) are dressed with their corresponding fractional PDF un¬ 
certainty (computed with each group’s dedicated set). This is 
a known fact and it has been extensively used in recent PDF 
benchmarks [48| . 

Given that correlations among the uncertainties are very 
small, as discussed above, and also that not every PDF set 
allows to take into account all the correlations, we define the 
combined PDF, a^ and mb uncertainty by simply summing the 
uncertainties in quadrature as 


^PDF + CKs+mf, 


\/+ {StY + 


( 11 ) 


Finally, we define the total theoretical uncertainty as the linear 
sum of the upper and lower variations for and (5pDF-i-aa+m6 
in a given flavour scheme. 

In table we report the NLO cross sections and their un¬ 
certainties at the 13-TeV LHC, for t-channel tH and tH pro¬ 
ductions separately, and for their sum tH + tH. Results are 
shown, using NNPDF2.3, in the 4F and 5F scheme for the 
static and dynamic scale choices in eqs. Q and 0- including 
the sources of uncertainty discussed above: scale uncertainty 
and combined PDF, as{mz) and mb one as well as the indi¬ 
vidual ones. The predictions in the combination of the 4F and 
5F schemes defined in eq. ([^ are presented in table The the¬ 
oretical uncertainty is dominated by the combined scale and 
flavour-scheme uncertainty 5 ^+fs over the PDF, Os and mb 
uncertainty 5pDF+as+mt ■ Figure [^summarises the NLO cross 
sections and the theoretical uncertainties for f-channel tH pro¬ 
duction, including the MSTW2008 and CTIO predictions. 

We conclude this section by commenting on two additional 
minor sources of uncertainty. The first one is related to the 
value of the Higgs and top quark masses. In tablej^we collect 
results for the t-channel NLO cross section (in the 5F scheme 
only) with parametric variations of 1 GeV in mn and mt- 
The variations have a modest impact on the total cross sec¬ 
tion, about 1 % only when both masses are varied in the same 
direction. From the combination of Tevatron and LHG exper¬ 
imental results 49 the top mass is currently known with a 


precision better than 1 GeV, while the combination of the lat¬ 
est ATLAS and CMS measurements of the Higgs mass gives 
a precision better than 0.5 GeV [^. We conclude that the 
impact of these uncertainties on the f-channel cross section at 
the LHC is below 1%. The last source of uncertainty we discuss 
is the Yukawa coupling of the bottom quark. We have checked 
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t-channel 


^(mo) 


[fb] 


a’'” 

'^PDF+cKs+mf, ^PDF 


M) 


[fb] 


^PDF+Qs+m^ ^PDF 


4F 


tH 

45.90(7) 

+3.6 

-6.3 

+2.3 

-2.3 

±0.9 

+0.6 

-0.9 

o q 

d ci 

+ 1 

46.67(8) 

+4.3 

-6.1 

+3.2 

-1.9 

±0.9 

+1.6 

-0.4 

+2.6 

-1.6 

tH 

23.92(3) 

+4.2 

-6.6 

+2.5 

-2.7 

±1.4 

+ 1.6 
-1.8 

+ 1.4 
-1.5 

24.47(5) 

+4.4 

-6.8 

+2.5 

-2.3 

±1.4 

+ 1.4 
-1.4 

+ 1.6 
-1.2 

tH + iH 

69.81(11) 

+3.2 

-6.6 

+2.8 

-2.5 

±0.9 

+ 1.6 
-1.7 

+2.1 

-1.6 

71.20(11) 

+4.3 

-6.5 

+3.0 

-2.4 

±0.9 

+2.0 

-1.1 

+2.0 

-1.9 


5F 


tH 

48.80(5) 

+ 7.1 
-1.7 

+2.8 

-2.3 

±1.0 

+ 1.7 
-1.1 

+2.0 

-1.8 

47.62(5) 

+7.4 

-2.2 

+3.0 

-2.4 

±1.0 

+1.6 

-0.8 

+2.4 

-2.0 

iH 

25.68(3) 

+6.8 

-2.0 

+3.4 

-2.9 

±1.4 

+ 1.9 
-1.5 

+2.5 

-2.0 

25.07(3) 

+7.4 

-2.1 

+3.2 

-2.9 

±1.4 

+1.7 

-1.8 

+2.4 

-1.8 

tH + iH 

74.80(9) 

+6.8 

-2.4 

+3.0 

-2.4 

±1.0 

+ 1.5 
-1.1 

+2.4 

-1.9 

72.79(7) 

+7.4 

-2.4 

+2.9 

-2.3 

±1.0 

+1.2 

-1.4 

+2.4 

-1.6 


Table 2. NLO cross sections and nncertainties for pp tHq^ iHq and {tHq + iHq) at the 13-TeV LHC. NNPDF2.3 PDFs 
have been nsed (NNPDF2.1 for mb uncertainty in 5F). The integration nncertainty in the last digit(s) (in parentheses) as well 
as the scale dependence and the combined PDF + + mb uncertainty in eq. (Ill (in %) are reported. The individual PDF, Qs 

and mb uncertainties are also presented as a reference. 


f-channel 

(mo) rTT] r% r% 

•^NLO ^Ai+FS ‘^PDF+as+mb 

(/^o) [fUl x% x% 

*^NLO L'-’J ^m+FS ^PDF+as+mb 

4F±5F tH 

iH 

tH + iH 

47.64(7) ±9.7 tli 

24.88(4) ±10.2 tli 

72.55(10) ±10.1 1^-4 

47.47(6) ±7.7 

24.86(3) ±8.3 12:3 

72.37(10) ±8.0 til 

Table 3. Same as table 

I 

but for the flavour-scheme combined results, according to eq. ^ 
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Fig. 4. Summary plot of the NLO cross sections with un¬ 
certainties for Higgs production associated with a single top 
quark, via a t-channel W boson, at the 13-TeV LHC. For the 
uncertainties, the inner ticks display the scale (plus combined 
flavour-scheme) dependence 5 ;j(+ps), while the outer ones in¬ 
clude (JpoF-t-Ca+mf 


that it is completely negligible, both in the 4F and 5F schemes, 
the impact of turning yt on/off at NLO being smaller than the 
numerical accuracy (0.1 —0.2%). Finally, we remind the reader 
that EW corrections for this process are presently unknown, 
and these could have an impact on the accuracy of the present 
predictions. 


[fb] 

172.3 

mt 

173.3 

174.3 

124.0 

75.54 (+1.0%) 

75.18 (±0.5%) 

74.99 (±0.3%) 

mn 125.0 

75.10 (±0.4%) 

74.80 

74.43 (-0.5%) 

126.0 

74.70 (-0.1%) 

74.16 (-0.8%) 

73.74 (-1.4%) 


Table 4. Higgs and top quark mass dependence of the NLO 
cross sections in the 5F scheme for pp —>■ tHq -|- iHq at the 
LHC with = 13 TeV. NNPDF2.3 PDFs have been used 
with po = {mu + The figures in parentheses are the % 

variations with respect to the reference cross section, computed 
with mu ~ 125.0 GeV and mt = 173.3 GeV. 


3.3 Distributions 

We now present a selection of kinematical distributions for 
the combined f-channel tH -|- iH production at the 13-TeV 
LHC, with NLO corrections and matching to a parton shower 
(NLO-I-PS). For the sake of brevity, we do not consider top 
and anti-top processes separately in this section, and will dub 
with t both the top quark and its antiquark. Our main interest 
here is to assess the precision of the predictions for f-channel 
production, therefore we do not specify any decay mode for the 
Higgs boson, i.e. we leave it stable in the simulation. On the 
other hand, we consider (leptonic) top decays, which allows us 
to compare the distributions of 6-jets coming from the hard 
scattering to the ones coming from the top quark. 

For the kinematical distributions, we use NNPDF 2.3 PDFs 
and the PythiaS parton shower. We have compared predic¬ 
tions obtained with the MSTW2008 and CTIO PDF sets and 
found no difference worth to report. We have also employed 
the HERWIG6 parton shower to verify that some important 
conclusions on the difference of the radiation pattern between 
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Fig. 5. Representative differential distributions for the Higgs boson and the top quark at NLO+PS accuracy in t-channel tH 
associated production at the 13-TeV LHC. The lower panels provide information on the differences between 4F and 5F schemes 
as well as the differential K factors in the two schemes. 


4F and 5F schemes were not dependent on shower programs. 
We estimate the scale dependence by varying and fx f inde¬ 
pendently by a factor two around the reference dynamic scale 
Ht!^ defined in eq. 0: which provides smaller scale depen¬ 
dence than the static choice for differential distributions, espe¬ 
cially for the high-pT region. 

We start by showing in fig. differential distributions for 
the Higgs boson and the top quark (before they decay). The 
first observation is that NLO distributions in the 4F and 5F 
schemes are in excellent agreement within their respective un¬ 
certainty associated to scale variation, i.e. within the 10% level. 
Interestingly, though, differential K factors (information in the 
insets below) are more pronounced for the 5F than for the 
4F scheme, the NLO results in the 5F scheme typically being 
out of the uncertainties as estimated from scale variation at 
LO. It should be noted that the LO process in the 5F scheme 


does not depend on the renormalisation scale, and therefore 
its smaller uncertainty (especially in the high-pr region) can 
be an artefact of the scheme. Results in the 5F tend to have a 
scale uncertainty that increases with pr much more than in the 
4F, but in most cases the differences are not striking. Slightly 
larger deviations between 4F and 5F appear only very close to 
the tH threshold, a region where we expect the 4F scheme to 
catch the underlying physics already at LO. 

In fig. § we present distributions for the two hardest jets 
which are not tagged as 6-jets. Jets and b-jets are defined in 
eqs. 0 and 1^. The contributions from the non-taggable for¬ 
ward 6-jets (2.5 < IpI < 4.5) are also denoted by shaded his¬ 
tograms as a reference. The jet with the highest transverse 
momentum (ji) tends to be produced in the forward region, 
very much like in single-top and VBF production. Most of the 
time this jet can be clearly associated to the light-quark current 
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Fig. 6. Same as 



shaded histograms as a reference. 


in the hard scattering. The very good agreement between 4F 
and 5F is manifest. This is expected as this observable should 
not be too sensitive on the details of heavy-quark current, as 
colour connections between the two currents are either vanish¬ 
ing or suppressed at the order in QCD we are working. On the 
other hand, sizeable differences arise for the second-hardest jet 
(^2), which shows a much steeper pr spectrum and tends to 
be produced centrally. The difference between predictions in 
the 4F and 5F schemes is often much larger than the scale un¬ 
certainty band (which is more pronounced in the 5F scheme 
in the bulk of the events). We will discuss further this feature 
when presenting jet multiplicities in the following. 

In fig. [7] we show the analogous distributions for the b- 
tagged jets. These are all the jets containing a 6-hadron and 
falling inside the acceptance of the tracking system, eq. ||^. 
We consider the two hardest b-jets {jb,i and 46,2) in the event 
regardless of their origin and, separately, we study the 6-jet 


coming from the top quark decay jb,t (tagged by using Monte 
Carlo information). The pr spectrum of jb,i has a rather long 
tail compared to jb ^2 and, at variance with light jets, all the 6- 
jets tend to be produced in the central region. Scale dependence 
at NLO is rather small in the 4F scheme, never reaching 10% 
and being typically around 5%. Differences between 4F and 5F 
predictions are visible, specially in the uncertainty band of jb ,2 
in the 5F scheme; this is of course expected, given that this 
observable is described only at LO accuracy in this scheme. 
Quite remarkably, however, these differences at NLO are often 
significantly less pronounced than in the case of light jets (spe¬ 
cially for the second jet), while naively one might expect the 
6-jet observables to be mostly affected by the flavour-scheme 
choice. On the other hand, at LO the inadequacy of the 5F 
scheme to describe 6-jets is evident. 

Comparing the transverse momentum of jb,t (first row, 
right plot in fig. to the corresponding spectra of jb,i and 
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Fig. 7. Same as fig.|^ but for the 6-tagged jets. On the right column the distributions for the &-jet coming from the top quark 
decay, selected by using Monte Carlo information, are shown. 


jb,2, it can be inferred that 6-jets from the top quark mostly 
contribute to the hardest 6-jet {jb,i) spectrum at low pr- On 
the other hand, as the pr tail falls much more rapidly for jb^t 
than for jb,i, gluon splitting in the hard scattering is the pre¬ 
dominant mechanism at high pr, and thus the main source of 
6-jets in this region. This observation also explains why the 
scale dependence in the 5F is small for low PT{jb,i), which is 
described at NLO accuracy, and increases sharply in the high- 
PT{jb,i) region, where the physics is dominated by the trans¬ 
verse dynamics of the g ^ bb splitting, which is described only 
at LO. 

We conclude this section by studying the jet multiplicities, 
which are sensitive to the flavour scheme as well as to the 
choice of the shower scale. As argued in [^, the dynamics of 
g —^ 66 splitting takes place at a scale which is typically lower 
than the hard scale of the process mt + mn or Ht, affecting 
the choice for the factorisation scale that one should use to 


describe t-channel production. An analogous argument could 
be made also for the shower scale choice 


15 


which in the 

MadGraph5_aMC@NLO matching procedure is chosen to be 
of the order of the partonic centre-of-mass energy in the Born 
process. In fig. we study the dependence of jet rates on 
the flavour scheme as well as on the shower scale, where two 
different choices of the shower scale are compared: one is the 
default value, and another is the default value divided by a 
factor of four. We can see that reducing the parton-shower 
scale has only a minor impact on the distributions, while a 
more interesting pattern arises from the choice of the flavour 
scheme. 


For the 6-tagged jets (right panel in fig. 1^, differences be¬ 
tween the two schemes are rather mild (~ 15% in the 2-jet bin 
and less for 0 and 1 jet) and always compatible within the scale 
uncertainty, which for the 2-jet bin is much larger in the 5F 
(the accuracy being only at LO). 
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t-channel tH at the LHC 13 
NLO+PythiaS 


da/dN(j) [fb/bin] 
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Fig. 8. Jet rates at NLO+PS accuracy in 4F and 5F schemes 
with different choices of the shower scales. 


Fig. 9. Jet rates only for the light jets both at fixed order 
and matched to a parton shower in 4F and 5F schemes with 
different choices of the shower scales. 


For non-fo-tagged jets (left panel in fig. |^, on the other 
hand, a higher jet multiplicity is clearly observed in the 4F 
scheme, which implies that harder QCD radiation is favoured 
in this scheme. Interestingly, the difference is visible already 
at the 1-jet bin, which is described at NLO accuracy at the 
matrix-element level. These differences cannot arise from the 
small component of forward, non-taggable heavy jets; on the 
contrary, they can be understood by considering jets that come 
from genuinely light QCD radiation. In fig. we show explicitly 
the multiplicity of light jets only (tagged by using Monte Carlo 
information), both at fixed order in QCD and at NLO matched 
to parton shower. Our first observation is that results in the 4F 
and 5F are almost identical at fixed LO (where only the zero 
and one jet bins are filled). The difference is therefore an effect 
of higher-order corrections, as it is confirmed by observing the 
fixed-NLO histograms. We recall that the fixed-order matrix 
element has a different colour structure in different schemes; 
in particular, the 4F at LO features a gluon in the initial state 
(compared to the 6-quark in the 5F) and an extra b in the 
final state. The radiation of extra light QCD partons from the 
g ^ bb splitting is therefore favoured in the 4F (e.g. an extra 
gluon can either attach to the initial-state gluon or to one of 
the 6’s, while in the 5F it can attach only to the initial-state 
6). This is indeed what we observe at fixed NLO. 

If the origin of the difference in the jet rates can be traced 
back to the difference between the LO 4F and 5F colour struc¬ 
tures, then one would also expect this difference to be miti¬ 
gated once higher-order corrections are included. To this aim, 
we have performed a fixed-order computation of the 2-jet bin in 
the 5F at NLO accuracy, i.e. calculated tHjj at NLO, within 
our simulation framework, finding indeed that the rate is sig¬ 
nificantly enhanced (by ~ 60%), lying much closer to the 4F 
result. A further hint that the scheme difference is indeed mit¬ 
igated at higher orders is given by the NLO-fPS results, which 
show that the 2-jet bin in the 4F is reduced by ~ 10% after the 
shower, while the corresponding 5F one is enhanced by ~ 30% 
over the fixed-order result. Finally, we have checked that the 
same results we have found here for single top plus Higgs, occur 


also in the case of single top production alone. In conclusion, 
our results suggest that the inclusion of the g —> 66 splitting 
in the matrix-element description at the lowest order, i.e. the 
4F scheme, allows a wider range of observables relevant for the 
analyses to be described more accurately. 


4 s-channel production 

Higgs-top quark associated production at hadron colliders can 
also be mediated by s-channel diagrams, see fig. Compared 
to f-channel production, the s-channel mechanism is naturally 
suppressed by the higher virtuality of the intermediate W bo¬ 
son and features a much smaller cross section at the LHC. In 
this section we calculate the NLO cross section, evaluating the 
corresponding uncertainties, and compare s-channel distribu¬ 
tions to those of f-channel production at NLO-I-PS level. 

At LO, s-channel production proceeds through qq annihi¬ 
lation into a virtual W boson, which can either emit a Higgs 
boson and then split to a tb final state, or first split to tb with 
the subsequent emission of a Higgs from the top quark. It turns 
out that in this case the interference between these two dia¬ 
grams is positive and its effect are much less relevant than in 
f-channel production [10| . At NLO, extra radiation can take 
place from either initial or final state, with no interference be¬ 
tween the two due to colour conservation. For the same reason, 
no interference between the s-channel and f-channel processes 
is present in the 5F scheme and the separation between chan¬ 
nels is still exact at NLO accuracy. In this production mode, 
bottom quarks are directly produced in the hard scattering 
via electroweak interaction and appear only in the final state. 
Thus, at variance with the f-channel and JF-associated pro¬ 
duction, the flavour scheme is not a key source of uncertainties 
for s-channel production. 

In the MadGraph5_aMC@NLO framework the code and 
the events for s-channel production at hadron colliders can be 
automatically generated by typing the following commands: 
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Fig. 10. Shape comparison between s- and t-channel distributions for the Higgs boson and the top quark at NLO+PS accuracy. 
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Fig. 11. In the top row: shape comparison between s- and f-channel distributions for jet rates (left), pr (centre) and p (right) 
spectra for the hardest jet at NLO-fPS accuracy. In the bottom row: corresponding plots for ^-tagged jets. 






































































12 


F. Demartin et al.: Higgs production in association with a single top quark at the LHC 


s-channel 

cnlo [fb] 

^11 

<^PDF 

5^” 

tH + iH 

2.812(3) 

+ 1.6 
-1.2 

+ 1.4 
-1.4 

+0.3 

-0.5 


Table 5. NLO total cross section for the processes pp —>■ tHb + 
iHb via an s-channel IT-boson exchange at the LHC {-y/s = 
13 TeV). NNPDF2.3 PDFs have been used. The integration 
uncertainty in the last digit (in parentheses), the fractional 
scale dependence and the PDF and Ofs uncertainties (in %) are 
also reported. 


(> import model loop_sm) 

> generate p p > w+ > h t b~ [QCD] 

> add process pp>w->ht~b [QCD] 

> output 

> launch 

In table[5]we show the total cross section at NLO. Reference 
values for the factorisation and renormalisation scales are set 
to po = Ht/2 = rriT 12 ■ Being a pure EW process at LO, s- 
channel production exhibits very low scale and uncertainties 
up to NLO. In the SM, the total rate amounts to about 3 fb, 
i.e. less than 5% of the t-channel cross section. 

In figs. |10|and [TT] we compare the shape of some distribu¬ 
tions between the s-channel and t-channel production modes at 
NLO-I-PS accuracy. We can see that most of the observables re¬ 
lated to s-channel events display a significantly different shape. 
Even though the total cross section in s-channel production is 
tiny and deviations from a f-channel-only simulation would 
probably fall inside the uncertainty band, the s-channel simu¬ 
lation can be included with little extra computing cost when 
precision is needed (it is also extremely fast at NLO). 


5 Higgs characterisation 


In this section we go beyond the SM and explore the sensitiv¬ 
ity of Higgs-single-top associated production to a Higgs boson 
coupling to the top quark that does not conserve CP. Several 
phenomenological studies on anomalous Higgs coupling deter¬ 
mination via Higgs-single-top associated production have ap¬ 
peared [^[^ |51||5^ . Current experimental constraints on the 
Higgs-boson couplings favour the SM, and in particular for the 
top quark the magnitude is consistent with the SM expecta¬ 
tions, even though an opposite sign with respect to the SM one 


is not yet completely excluded 57 58 


Moreover, although the scenario of a pseudoscalar Higgs is 
disfavoured |59| 60 , no stringent constraint has been put on 


a CP-violating Htt coupling. In fact, even if current results 
are fully compatible with the SM hypothesis, some analyses on 
public LHC data seem to favour a non-zero phase in the top 
quark Yukawa interaction [61}j64| . 

In this work we consider the (simplified) case of a spin-0 
particle with a general CP-violating Yukawa interaction with 
the top quark, which couples both to scalar and pseudoscalar 
fermionic densities. On the other hand, we assume the inter¬ 
action with the W bosons to be the SM one. We note that 
this assumption does not correspond to a typical realisation of 
CP-violation in a two-Higgs-doublet model (2HDM) where the 
mass eigenstates are CP-mixed states and their coupling to the 
vector bosons is reduced. Our setup, however, corresponds to 



Fig. 12. NLO cross sections (with scale uncertainties) for ttXo 
and f-channel tXo productions at the 13-TeV LHC as a func¬ 
tion of the CP-mixing angle a, where khu and are set to 
reproduce the SM GF cross section for every value of a. 


considering the effective SM Lagrangian and to including the 
operator 


C = ^('('V) Qh^tR + h.c. 


( 12 ) 


with ct complex. The implementation we use is based on the 
effective field theory framework presented in refs. |24|j26] and 
employs the HC_NLO_XO model 


27 


^^The effective Lagrangian 
for the Higgs-top quark interaction (|12 1 below the EWSB scale 
leads to (see eq. (2.2) in ref. [24|) 


To = —1pt{CaKHttgHtt + iScKAttgAtt 75)V’t , 


(13) 


where Xq labels a generic spin-0 particle with CP-violating 
couplings, Ca = cos a and Sa = sin a are related to the CP- 
mixing phase a, K,Htt,Att are real dimensionless rescaling pa¬ 
rameters, and gntt = gAtt = mt/v (= ytl\/2), with v ~ 246 GeV. 
While redundant (only two independent real quantities are 
needed to parametrise the most general CP-violating interac¬ 
tion with the top quark at dimension four), this parametrisa- 
tion has the practical advantage of easily interpolating between 
the CP-even (cq = l,Sa = 0) and CP-odd (cq = 0, Sa = 1) 
couplings, as well as to easily recover the SM case by setting 
Ca = 1 , = 1. 

The nature of the top quark Yukawa coupling directly af¬ 
fects the loop-induced Higgs coupling to gluons (together with 
an effect on the couplings to 77 and Z'^, which are also modi¬ 
fied but not considered here) 

fyg I / ^a,fj,iy 

^0 ^ [Ca^rfgg^?/fgg ^ 

+ SaACAggSAgg , (14) 

^ For the code and event generation, one can simply issue the 
command ‘import model HC_NL0_X0’ and replace ‘h’ by ‘xO’ in 
MadGraph5_aMC@NLO. 
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Fig. 13. Differential distributions for the Higgs boson and the top quark at NLO+PS accuracy in t-channel 
production at the 13-TeV LHC, with different values of the CP-mixing angles, where Khu and kah are set 
reproduce the SM GF cross section for every value of a. 


tH associated 
in eq. (161 to 




Fig. 14. Shape comparison among different values of the CP-mixing angles, where Khh and Kah are set in eq. | |16[ ) to reproduce 
the SM GF cross section for every value of a. Pseudorapidity separation between the Higgs and the top quark (left) and opening 
angle between the hardest jet and the lepton from the top quark in the lab frame (right). 


where Qngg = —ols!(Sttv) and gAgg = as/{2'iTv). In the parametri- 
sation given above, the strength of the coupling between Higgs 
and gluons can be rescaled independently of the top quark 
Yukawa coupling. Assuming that the the top quark dominates 
the gluon-fusion (GF) process at the LHG energies, then K.Hgg —>■ 

, K,Agg “> ^Att ■ lu SO doing, the ratio between the actual 
cross section for GF at NLO QCD and the corresponding SM 
prediction can be written as 


„93^Xo 

"nlo 

gg—^H 

^NLO,SM 


2 2 
Cq: 


+ S( 


: I f^Att - I 

^ QHgg ' 


(15) 


because there is no interference between the scalar and pseu¬ 
doscalar components in the amplitudes for Higgs plus up to 
three external partons, see e.g., [^. In particular, if the rescal¬ 
ing parameters are set to 


f^Htt 1 , I^Att I QHgg/QAgg | 2/3 , 


(16) 


the SM GF cross section is reproduced for every value of the 
GP-mixing phase a. Given that current measurements are com¬ 
patible with the expected SM GF production rate, one can 
consider the simplified scenario where the condition in eq. (161 
is imposed and the CP-mixing phase a is basically left uncon¬ 
strained by current data. 


Figure shows the total cross section for t-channel tXo 
production as a function of the CP-mixing angle a. We also 
show the ttXo cross section, which is not only another pro¬ 
cess sensitive to the modifications of the top quark Yukawa 
coupling in eq. (131, but also a background to t-channel pro¬ 
duction. The uncertainty band represents the envelope defined 
in sect. |3.2| i.e. the combined scale and flavour-scheme depen¬ 
dence. The ttXo uncertainty band represents the scale depen¬ 
dence only, when the scale is va ried by a factor two around 


MO = mT{t) mriXo) 


26 


The first important observation is that while the GF and 
ttH cross sections are degenerate under yt —> —yt (depend¬ 
ing quadratically from the top quark Yukawa coupling), in 
t-channel production this degeneracy is clearly lifted by the 
interference between diagrams where the Higgs couples to the 
top quark and to the W boson. In it was shown that the 
t-channel cross section is enhanced by more than one order of 
magnitude when the strength of the top Yukawa coupling is 
changed in sign with respect to the SM value. Here we can see 
how the same enhancement can take place also in the pres¬ 
ence a continuous rotation in the scalar-pseudoscalar plane. 
While not affecting GF (by construction), such a rotation has 
an impact also on the ttXo rate, which is in general lower 
for a pseudoscalar or CP-mixed state [26|. t-channel produc- 
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tion lifts another degeneracy present in GF and ttXo, namely 
a —>■ TT — a. Given the partial compensation between the t- 
channel and tiXo cross sections at different values of q, an 
analysis which could well separate between the two production 
mechanisms would be needed to put stringent constraints on a 
CP-violating Higgs coupling to the top quark. 

We remind that the enhancement of the t-channel cross sec¬ 
tion takes place mostly at threshold, as one can clearly see in 
the left plot of fig. |13| This means that one should not be con¬ 
cerned by violations of perturbative unitarity at the LHC, as 
they do not appear for partonic centre-of-mass energies lower 
than ~ 10 TeV . In fig. we also show the transverse mo¬ 
mentum distributions for the Higgs and the top quark. The 
distributions are well behaved in this case too, not displaying 
any strong trend in their high-pT tails, i.e anything that could 
suggest a unitarity violating behaviour. 

Finally, in fig. we plot the pseudorapidity separation 
between the Higgs and the top quark (left) and the opening 
angle between the hardest jet and the lepton from the top 
quark in the lab frame (right), showing that these variables 
have a discriminating power on a. For this last observable, the 
lepton is required to satisfy the following selection criteria 

PT{e) > 20 GeV , \rj{£)\ < 2.5 . (17) 

6 Summary 

In this work we have studied the production of a Higgs boson 
in association with a single top quark at the LHC. Our aim 
has been to carefully consider the effects of NLO corrections in 
QCD on total cross sections and differential distributions for 
t- and s-channel production. We have scrutinised a wide range 
of theoretical systematic uncertainties and in particular those 
arising from the choice of the heavy-quark scheme, 4-flavour or 
5-flavour. We have found that at the level of total cross sections 
a comfortable consistency between the two schemes exists when 
physically motivated choices for the renormalisation and fac¬ 
torisation scales are made, with similar resulting uncertainties. 
For differential distributions, on the other hand, the situation is 
slightly more involved. While sizeable differences between the 
two schemes arise at LO, they are considerably milder at NLO 
and NLO-I-PS, in line with expectations. In this case, we have 
shown that the 4F and 5F schemes provide fully consistent and 
similarly precise predictions for distributions such as those of 
the Higgs boson, the top quark, and the forward jet. On the 
other hand, the 4-flavour scheme is in general able to provide 
accurate predictions for a wider set of observables, including 
those of the spectator fe-quark and extra jets. In addition to 
t-channel production in the SM, we have also briefly presented 
the results for the subdominant s-channel production, high¬ 
lighting the differences in the most important distributions 
with respect to the corresponding ones of t-channel produc¬ 
tion. Finally, we have provided results (total cross sections as 
well as a few representative distributions) for the case where an 
explicit CP violation is present in the coupling between the top 
quark and the Higgs boson, making it clear that in this case 
Higgs associated production with a single top could provide 
complementary and very valuable information to that of tiH 
production. We conclude by stressing that all results presented 
here have been obtained by employing the publicly available 
MadGraph5_aMC@NLO framework and therefore they can 
be easily reproduced (and possibly extended) by generating 


the corresponding event samples to be used in fully-fledged 
experimental analyses. 
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